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ABSTRACT 

We study numerically the thermal emission of e + e~ pairs from a bare strange 
star heated by energy input onto its surface; heating starts at some moment, and 
is steady afterwards. The thermal luminosity in e + e~ pairs increases to some 
constant value. The rise time and the steady thermal luminosity are evaluated. 
Both normal and colour superconducting states of strange quark matter are con- 
sidered. The results are used to test the magnetar model of soft gamma-ray 
repeaters where the bursting activity is explained by fast decay of superstrong 
magnetic fields and heating of the strange star surface. It is shown that the rise 
times observed in typical bursts may be explained in this model only if strange 
quark matter is a superconductor with an energy gap of more that 1 MeV. 

1. Introduction 

Strange stars are astronomical compact objects which are entirely made of deconfined 
quarks (for a review, see Glendenning 1996; Weber 1999). The possible existence of strange 
stars is a direct consequence of the conjecture that strange quark matter (SQM) may be 
the absolute ground state of the strong interaction, i.e., absolutely stable with respect to 
56 Fe (Bodmer 1971; Witten 1984). SQM with a density of ~ 5 x 10 14 g cm -3 might exist up 
to the surface of a strange star. Recently, the thermal emission from bare SQM surfaces of 
strange stars was considered (Usov 1998, 2001a). It was shown that the surface emissivity 
of SQM in both equilibrium photons and e + e~ pairs created by the Coulomb barrier at 
the SQM surface is > 10% of the black body surface emissivity at the surface temperature 
T s > 1.5 x 10 9 K. Below this temperature, T s < 1.5 x 10 9 K, the SQM surface emissivity 
decreases rapidly with decrease of T s . 
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At the moment of formation of a strange star the surface temperature may be as high 
as ~ 10 11 K (e.g., Haensel, Paczyhski, & Amsterdamski 1991). Since SQM at the surface 
of a bare strange star is bound via strong interaction rather than gravity such a star can 
radiate at the luminosity greatly exceeding the Eddington limit, up to ~ 10 52 ergs s _1 
at T s ~ 10 11 K. (Alcock, Farhi, & Olinto 1986; Chmaj, Haensel, & Slomihski 1991; Usov 
1998, 2001a). A young strange star cools rapidly, and within about a month after its 
formation the surface temperature is less than 2.5 x 10 8 K (e.g., Pizzochero 1991). In this 
case, the thermal luminosity from the stellar surface in both equilibrium photons and e + e~ 
pairs is negligibly small, L t h < 10 25 ergs s _1 . A bare strange star with such a low surface 
temperature may be a strong source of radiation only if its surface is reheated. Recently, 
response of a bare strange star to accretion of a massive comet-like object with the mass 
AM ~ 10 25 g onto the stellar surface was considered (Usov 2001b). It was shown that the 
light curves of the two giant bursts observed from the soft 7-ray repeaters SGR 0526-66 
and SGR 1900+14 may be easily explained in this model. 

In this paper we report on our numerical simulations of the response of a bare strange 
star to energy input onto its surface. We consider a wide range of the rate of the energy 
input. Both normal and colour superconducting SQM are discussed. 

2. Formulation of the problem 

The model to be studied is the following. The energy input onto the surface of a bare 
strange star starts at the moment t — 0, and it is spherical and steady at t > 0. Since in our 
simulations the surface temperature is not higher than a few x 10 9 K, e + e~ pairs created 
by the Coulomb barrier are the main component of the thermal emission from the stellar 
surface (Usov 2001a). We assume that the process of the energy input has no effect on the 
outflow of both created pairs and photons which form due to annihilation of some of these 
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pairs (cf. Usov 2001b). 

The equation of heat transfer that describes the temperature distribution at the surface 
layers of a strange star is 

c <*L = JL( K —)-£ (i) 

dt dx \ dx ) 

where C is the specific heat for SQM per unit volume, K is the thermal conductivity, and 
e v is the neutrino emissivity. 

The heat flux due to thermal conductivity is 



q = -KdT/dx . (2) 

At the stellar surface, x — 0, the heat flux directed into the strange star is equal to (Usov 
2001b) 

q = L- mput /(4nR 2 )-e ± f ± , (3) 

where L mpVit is the rate of the energy input onto the stellar surface, R ~ 10 6 cm is the 
radius of the star, e±f± is the energy flux in e + e~ pairs emitted from the SQM surface, 
e± ~ m e c 2 + kT s is the mean energy of created pairs, 



f ± ~ 10 39 - 2 T s 3 9 exp (-±j£) J(C) crn^s" 1 (4) 

is the flux of pairs from the unit SQM surface, 



1 C 3 In (1 + 2C- 1 ) vr 5 C 4 ( . 

lU 3 (1 + 0.074C) 3 6 (13.9 + C) 4 ' 
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(~(2x 10 10 K)/T s , and T Si g is the surface temperature in units of 10 9 K. 

Eqs. (2) — (5) give a boundary condition on dT/dx at the stellar surface. We assume 
that at the initial moment, t — 0, the temperature in the surface layers is constant, 
T = 10 8 K. 

It has been suggested (Bailin & Love 1979, 1984) that the quarks may eventually form 
Cooper pairs. Recently, superconductivity of SQM was considered in detail (for a review, 
see Rajagopal, & Wilczek 2000; Alford, Bowers, & Rajagopal 2001), and it was shown 
that SQM is plausibly a colour superconductor if its temperature not too high. Below, we 
consider both normal and superconducting SQM. 

2.1. The Normal State of SQM 

For non-superconducting SQM, the contribution of the quarks to both the specific 
heat and the thermal conductivity prevails over the contributions of the electrons, photons 
and gluons. In this case we have (Iwamoto 1982; Heiselberg & Pethick 1993; Benvenuto & 
Althaus 1996) 

C ~ C q ~ 2.5 x 10 20 (n fe /n ) 2/3 T 9 ergs cm" 3 K" 1 (6) 

K ~ K q ~ 6 x 10 20 a-\n b /n ) 2/;i ergs cm" 1 s" 1 K" 1 (7) 

e u ~ 2.2 x 10 26 a c F e 1/3 (n 6 /n )T 9 6 ergs cm" 3 s" 1 (8) 

where uq ~ 1.7 x 10 38 cm" 3 is normal nuclear matter density, n\, is the baryon number 
density of SQM, a c = g 2 '/4ii is the QCD fine structure constant, g is the quark-gluon 
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coupling constant, Y e = n e /rib is the number of electrons per baryon, and Tg is the 
temperature in units of 10 9 K. 

2.2. The Colour Superconducting State of SQM 

SQM may be a colour superconductor if its temperature is below some critical 
value. In the classic model of Bardeen, Cooper, and Schrieffer the critical temperature is 
T c ~ 0.57A /A: B , where A is the energy gap at zero temperature and k B is the Boltzmann 
constant (e.g., Carter & Reddy 2000). The value of A is very uncertain and lies in the 
range from — 0.1 — 1 MeV (Bailin & Love 1984) to ~ 50 - 10 2 MeV (Alford, Rajagopal, & 
Wilczek 1998; Alford, Berges, & Rajagopal 1999; Pisarski k Rischke 2000). 

We use the following interpolation formula for the temperature dependence of the 



energy gap(e.g., Carter & Reddy 2000): 



A(T) = A c 



1 T, 



rji \ 2 



1/2 



(9) 



Superconductivity modifies the properties of SQM significantly. The specific heat of 
superconducting SQM increases discontinuously as the temperature falls below the critical 
temperature, and then decreases exponentially at lower temperatures. For the quark specific 
heat at T < T c we adopt the theoretical results of Miihlschlegel (1959) for superconducting 
nucleons, i.e., (Horvath, Benvenuto, & Vucetich 1991; Blaschke, Klahn, & Voskresensky 
2000) 

C q = 3.2C q (T c /T) exp (-A /T) 
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x [2.5 - l.7T/T c + 3.6(T/T C ) 2 ] , (10) 

here and below tilde signifies that this value relates to superconducting SQM. Even 
at T C T c the suppression of the specific heat of SQM is never complete because the 
electrons remain unpaired. The specific heat of superconducting SQM which is used in our 
simulations is C ~ C q + C e , where 

C e ~ 5.7 x 10 19 F e 2/3 K/n ) 2/3 T 9 ergs cm" 3 K" 1 , (11) 

is the specific heat of the electron subsystem of SQM (Blaschke, Grigorian, & Voskresensky 
2001). 

At T < T c both the thermal conductivity of SQM and its neutrino emissivity are 
suppressed by a factor of exp (—A/k B T). In our simulations, we adopt 

K = K q exp[-A(T)/k B T] + K e , (12) 

e„ = e„exp[-A(r)/A; B r], (13) 

where K q and e u are given by equations (7) and (8), respectively, and 

K e ~ 5.5 x 10 23 y e (V^o)r 9 _1 ergs cm" 1 s" 1 K" 1 (14) 

is the thermal conductivity of the electrons (e.g., Blaschke et al. 2001) 



3. Results of numerical simulations 

The set of equations (1) — (5) was solved numerically for both normal (Ao = 0) and 
superconducting (A > 0) states of SQM. We assumed the typical values of a c = 0.1, 
rib — 2no, and Y e = 10~ 4 . Figure 1 shows a typical temporal behaviour of the strange star 
luminosity, L± = 4:irR 2 e±f±, in e + e~ pairs. From this Figure, we can see that L± increases 
eventually to its maximum value L± ax . The rate of this increase may be characterized by 
the rise time r± that is determined as a time interval from the initial moment, t — 0, to 
the moment when L± is equal to (l/2)L™ ax . The results of our simulations are presented 
in Tables 1 and 2. There is a critical value, L cr , of the input luminosity at which the 
dependence of L™ ax on Li nput changes qualitatively (see Table 1). L CT is ~ 10 40 ergs s _1 for 
normal SQM and ~3x 10 38 ergs s _1 for superconducting SQM with Ao > 0.3 MeV. At 
-^mput <) a few x L CT , L™ ax is about L input while when -L input is a few times below than L cr 
the thermal emission of e + e" pairs from the stellar surface is negligible. 

In our simulations, the rise time r± varies in a very wide range from ~ 10 6 s at 
^input < 10 40 ergs s- 1 and A < 0.1 MeV to ~ 10~ 8 s at L input ~ 10 45 ergs s" 1 and 
A > 1 MeV (see Table 2). 

At A > 1 MeV the thermal emission of e + e~ pairs does not depend on A because in 
this case both the specific heat of the quarks and their thermal conductivity are strongly 
suppressed, and the heat transport is mostly determined by the electron subsystem of SQM. 

The rise time of the luminosity in neutrinos is many orders of magnitude larger that 
r±, especially when r± is small. In our model the neutrino luminosity may increase up to 
Anput — L™ ax when t goes to infinity. This is because all energy which is delivered onto the 
stellar surface is radiated either from the surface by e + e~ pairs or from the stellar interior 
by neutrinos. 
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4. Discussion 

Since bare strange stars can radiate at the luminosities greatly exceeding the Eddington 
limit (Alcock et al. 1986; Chmaj et al. 1991; Usov 1998, 2001a), these stars are reasonable 
candidates for soft 7-ray repeaters (SGRs) that are the sources of brief (~ 10 -2 — 10 2 s) 
bursts with Super-Eddington luminosities, up to ~ 10 43 — 10 45 ergs s -1 . The bursting 
activity of a SGR may be explained by fast heating of the bare, rather cold (T s < 10 8 K) 
surface of a strange star up to the temperature of ~ (1 — 2) x 10 9 K and its subsequent 
thermal emission (Usov 2001a). The heating mechanism may be either fast decay of 
superstrong (~ 10 14 — 10 15 G) magnetic fields (Usov 1984; Duncan & Thompson 1992; 
Paczyhski 1992; Thompson & Duncan 1995; Cheng & Dai 1998; Heyl & Kulkarni 1998) or 
impacts of comet-like objects onto the stellar surface (Harwitt & Salpeter 1973; Newman & 
Cox 1980; Zhang, Xu, & Qiao 2000; Usov 2001b). The magnetar model of SGRs which is 
based on the first mechanism is most popular now. In this model, the magnetic energy of 
a strongly magnetized strange star (magnetar) may be released from time to time due to 
MHD instabilities. A violent release of energy inside a magnetar excites surface oscillations 
(e.g., Thompson & Duncan 1995). In turn, this shaking may generate strong electric fields in 
the magnetar magnetosphere which accelerate particles to high energies. These high energy 
particles bombard the surface of the strange star and heat it. At the input luminosities 
of ~ 10 41 — 10 45 ergs s -1 , which are typical for SGRs, the efficiency of reradiation of the 
partical energy by the stellar surface to e + e~ pairs is very high, L±/ L- mput ~ 1, especially 
for superconducting SQM (see Table 1). Since for SGRs the burst luminosities are at least 
a few orders of magnitude higher than L* ~ 4irm c c 3 R/(T T ~ 10 36 ergs s -1 , the outflowing 
e + e~ pairs mostly annihilate in the vicinity of the strange star (e.g., Beloborodov 1999). 
Therefore, at L mput ~ 10 41 — 10 45 ergs s _1 far from the star the luminosity in X-ray and 
7-ray photons practically coincides with L input , L 7 ~ L± — L* ~ L± ~ L input . 
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Two giant bursts were observed on 5 March 1979 and 27 August 1998 from 
SGR 0526 — 66 and SGR 1900 + 14, respectively. The peak luminosities of these bursts 
were ~ 10 45 ergs s _1 (Fenimore, Klebesadel, & Laros 1996; Hurley et al. 1999). In this 
case, from Table 2 the rise time expected in our model is < 10~ 3 s that is consistent with 
available data on the two giant bursts (Mazets et al. 1999). This is valid irrespective of that 
SQM is a colour superconductor or not. For typical bursts of SGRs the luminosities are 
~ 10 41 — 10 42 ergs s _1 (e.g., Kouveliotou 1995), and the observed rise times (~ lO^ 1 — 10~ 3 s) 
may be explained in our model only if SQM is a superconductor with the energy gap 
A > 1 MeV (see Table 2). 

This work was supported by the Israel Science Foundation of the Israel Academy of 
Sciences and Humanities. 
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Table 1: The maximum luminosity L™ ax in e + e pairs 
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10 44 


9.999 x 


10 44 
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Note. - - Units of L™ ax are ergs s 1 . In all cases when in this Table L± ax is equal to Li np ut we have 

(-L'input — L± )/iinput < 10 
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Table 2: The rise time r± for the thermal emission of e + e pairs 
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Note. — Units of t± are seconds. 
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Fig. 1. — The thermal luminosity in e + e pairs as a function of time for Li nput = 10 41 ergs s x 
and A = 1 MeV. 



